Dynamics in shear flow studied by X-ray Photon Correlation Spectroscopy 
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X-ray photon correlation spectroscopy was used to measure the diffusive dynamics of colloidal 
particles in a shear flow. The results presented here show how the intensity autocorrelation functions 
measure both the diffusive dynamics of the particles and their flow-induced, convective motion. 
However, in the limit of low flow/shear rates, it is possible to obtain the diffusive component of 
the dynamics, which makes the method suitable for the study of the dynamical properties of a 
large class of complex soft-matter and biological fluids. An important benefit of this experimental 
strategy over more traditional X-ray methods is the minimization of X-ray induced beam damage. 
While the method can be applied also for photon correlation spectroscopy in the visible domain, our 
analysis shows that the experimental conditions under which it is possible to measure the diffusive 
dynamics are easier to achieve at higher q values (with X-rays). 



I. INTRODUCTION 



II. DESCRIPTION OF THE EXPERIMENT 
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X-ray photon correlation spectroscopy (XPCS) offers 
an unique way to perform direct measurements on the 
mesoscale dynamics in a large class of complex soft- 
matter systems (see for example Refs. [2, [H, H, 0, H, IH ) . 
Based on the same principle as Dynamic Light Scattering 
(DLS) [7], XPCS is however not limited to non-turbid 
samples, multiple scattering can usually be neglected, 
and can access larger values of the momentum transfer q, 
probing smaller spacial distances. Two major difficulties 
associated with XPCS are the small scattering strength 
of most samples and their degradation induced by the X- 
rays. Using a high-brilliance 3 rd generation synchrotron 
source helps addressing the former issue but increases at 
the same time the importance of the latter. In the case 
of fluid samples, performing XPCS experiments under 
continuous flow can limit the effects of beam damage. In 
addition, fluid mixers, in which the temporal dependence 
of a process taking place in a flow device is mapped into 
a spatial dependence along the flow channel, can be used 
to perform time-resolved studies [1, Q . 



In the experiments reported here, XPCS was used to 
study the Brownian dynamics of a colloidal suspension 
of (weakly to non-interacting) hard spheres under shear 
flow. In general, the correlation functions measured by 
XPCS are affected by both, the flow-induced dynamics 
as well as the diffusive dynamics of the particles. While 
measuring the flow velocity profile can be a legitimate 
goal in many applications (see for example Ref. [lCj), the 
aim of the present experiment was to describe the diffu- 
sive component of the dynamics. The results presented 
here show that for a laminar flow, and in the limits of low 
shear rates, it is possible to decouple the components due 
to diffusive and convective motion of the particles. 



A. X-ray photon correlation spectroscopy 

The experiment was performed at the Troika beamline 
(ID10A) of the European Synchrotron Radiation Facility 
(ESRF) in Grenoble, France. A single bounce Si(lll) 
crystal monochromator was used to select 8keV X-rays, 
corresponding to a wavelength of A = 1.55 A. The rel- 
ative bandwidth was AA/A ~ 10~ 4 . A Si mirror down- 
stream of the monochromator was suppressing higher or- 
der frequencies. The source size (FWHM) was approx- 
imately 30 x 900 fim 2 (vxh) and the source-sample dis- 
tance was 46 m. In the vertical direction, the X-ray beam 
was focused by a beryllium compound refractive lens in 
order to enhance the flux. A transversely partially coher- 
ent beam was defined by slit blades with highly polished 
cylindrical edges. The slit size was 10 /im in the vertical 
and horizontal directions. A set of guard slits, placed 
just upstream of the sample was used to minimize the 
parasitic X-ray scattering from the beam-defining slits 
(see Fig.QJi). Under these conditions, the incident X-ray 
flux was ~ 5 • 10 9 photons/s. The scattering from the 
colloidal suspension was recorded with a Cyberstar scin- 
tillator point detector placed at a distance of 2.30 m from 
the sample. The detection area was typically limited to 
100 x 100 /im 2 by precision slits placed in front of the 
detector, corresponding to a few speckles. 

As it will be seen, the dynamics measured in a flow- 
ing sample is not isotropic. In the experiments described 
here, measurements were performed at momentum trans- 
fers covering a range of q-a from 1.5 to 5.7 (a is the radius 
of the particles) in "transverse flow" and "longitudinal 
flow" geometries (Fig. Q]i) . 

Normalized intensity fluctuation autocorrelation func- 
tions, 
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FIG. 2: Static SAXS scattering from the latex particles and 
fitted form factor for the solid spherical particles (solid line). 
The straight line with slope -4, representing Porod's law is 
a guide for the eyes. Inset: estimate of the structure factor 
S(q) obtained by dividing the total scattered intensity to the 
fitted value for the form factor. 



FIG. 1: (a) XPCS experimental setup for the flow geometry 
considered here. Longitudinal flow scans were performed by 
moving the detector in the vertical z direction, while trans- 
verse flow scans were performed by keeping the detector in 
the horizontal scattering plane and moving it along the x di- 
rection, (b) Schematics of the flow channel with quadratic 
section (edge length 2R). The velocity profile has a parabolic 
shape; h is the transverse size of the X-ray beam. 



low (Re < 10 3 ), the flow is in a laminar regime, which 
was well the case in all the measurements performed here, 
as can be seen in table |H Assuming a no-slip boundary 
condition at the wall of the flow cell, the velocity profile 
can therefore be described by 



(2) 
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were obtained using a "Flex" hardware correlator man- 
ufactured by Correlator.com that was connected to the 
output of the X-ray detector. 



Here (see Fig. [T|b), r is the distance from the center of 
the capillary and vq is the average flow velocity (defined 
using the average flow rate Q — AR 2 vq). The mean shear 
rate is defined as 7 = =^-. 



B. Flowcell 

The in-house made flowcell consisted of an aluminum 
sheet through which the flow channel (see sketch in 
Fig. |T|b) was cut by high-precision machining. Polymer- 
coated kapton foils were fusion-bonded to the face of the 
aluminum and held in place by a "sandwich-type" pres- 
sure applying system. A high-precision syringe pump 
obtained from Harvard Apparatus was used to push the 
fluid through the flow device. The width of the chan- 
nel (2R = 1000 /im) was much larger than the one of 
the beam (h — 10 /mi). Therefore, the scattering vol- 
ume could be approximated by a one-dimensional object. 
This made it possible to apply the model of a Poiseuillc 
flow - a Newtonian fluid performing a shear flow in a 
cylindrically-shaped flow device. 

The laminar character of the flow is determined by the 
Reynolds number which measures the ratio of the inertial 
to viscous forces, and is defined as Re = Halgg [j| Here 
2R is the characteristic length of the system (confer to 
Fig. [TJs), vq is the flow velocity of the fluid and 77 the 
dynamic viscosity. As long as the Reynolds numbers are 



C. Sample 

Charge stabilized latex spheres with a nominal radius 
of a=110nm dispersed in water were obtained from Duke 
Scientific Co. at a concentration of 10% per weight. Glyc- 
erol with a purity of >99% was obtained from Sigma- 
Aldrich Co. The spheres suspension was mixed with the 
glycerol and placed in a dessicator which was evacuated 
on top of a heatable magnetic stirrer. The temperature 
was set to approximately 75° C which facilitated water 
evaporation and decreased the viscosity. These condi- 
tions were applied for about two days, leaving the sam- 
ples not only essentially water-free but also de-gased, 
avoiding thus an important reason for bubble formation. 
The amount of suspension and glycerol was chosen such 
that the final volume fraction was $ = 0.1. 

The small angle X-ray scattering (SAXS) pattern of 
the sample is shown in Fig. O The time-averaged data 
was fitted, for values of qa > 5.0, to the form factor of 
spherical particles with the size distribution given by a 
Schultz distribution function, yielding a radius of a — 
1070 A and a size polydispersity < 7%. Also shown is 
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the static structure factor S(q), exhibiting the weakness 
of interactions (S(q) « 1 for most of the g-range). 



III. THEORY: XPCS IN A LAMINAR FLOW 

The central quantity describing the dynamics of the 
colloidal suspension is the dynamic structure factor, or 
intermediate scattering function (ISF), which in the case 
of non-interacting (statistically independent) particles is 
given by Q 



j v 

9l (q,t) cx (E* k (0)Ek(t)exp [-iq- (r\.(0) - r' k (t))] ) . 
k=l 

(3) 

Here q is the scattering vector, Ek{t) is the amplitude 
of X-rays scattered by particle k at time t and r'k(t) is 
the particle position in the presence of flow at time t. 
The summation is over all particles in the scattering vol- 
ume for a particular realization of the experiment and the 
ensemble average is performed over many different real- 
izations. In order to differentiate between the diffusive 
and the flow-induced particle displacements, it is useful 
to write 

exp \-iq- (r' k (0)-r'k(t)) 

exp [-iq- (ffc(O) - r fe (t))] ■ exp(icf • v t) ■ exp(icf • Svt). 

(4) 

Here fk(t) is the particle position due to only to the dif- 
fusive motion, vq is the average flow velocity throughout 
the sample and Sv is the flow velocity difference between 
positions fk(t) and ffe(O). 

By analyzing Eqs. [3] and 2J it can be concluded that 
the ISF depends on three main factors: 

i) the particle diffusion, with a characteristic time scale 
to = (Dq 2 )^ 1 (D is the diffusion coefficient, see below) 

ii) the amplitude factor, which measures the intensity 
scattered by each particle, Ik = (E^Q) ■ Ek(0)), but also 
the transit time of the particles through the scattering 
volume, tt — h/vo (h is the transverse beam size), 

iii) shear-induced effects which are described, as it will 
be shown, by a time scale t$ oc (q'jRcoscf))^ 1 , which 
is inversely proportional to the magnitude of the local 
velocity gradient j. Here <j> is the angle between the 
scattering vector and the local velocity v. 

Note that the importance of the fourth factor, namely 
the average velocity-dependent exp [iq - vot], was not em- 
phasized because it can only be detected in heterodyne 
correlation functions [Tlj . A homodyne PCS experiment 
(like the one reported here) detecting the square modu- 
lus of the ISF, |<7i(q, t)\ 2 can measure velocity gradients 
but not absolute velocities 1131. 

It can be demonstrated [13| that, if the three time 
scales described above are sufficiently different from each 
other, the ISF can be factorized in independent contribu- 
tions from the three relaxation mechanisms. The physical 



argument is that if one of the contributions (e.g. diffu- 
sion) to the ISF varies on a much faster time scale than 
the others, then the particles explore the phase space as- 
sociated with this particular relaxation mechanism while 
the other factors are (approximately) constant. In such 
a situation, the ISF under shear flow can be written as 

gi (q,t) oc exp(iqv t) {E*(0)E(t)) x 

( J>xp [-iq- (r fc (0) - r k (t))] \ x 
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exp [iq - Sv(k, l)t] 
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(5) 



In writing Eq. [5J we have also assumed that all parti- 
cles are identical scatterers and that, due to the fact that 
particles are indistinguishable and that the system is sta- 
tionary, the shear induced effects can be calculated also 
by taking into account the velocity differences between 
all pair of particles for a particular realization of the sta- 
tistical ensemble. 

The physical interpretation of the shear-induced effects 
is that the frequencies of X-rays scattered by particles 
which are moving with different (average) velocities are 
Doppler-shifted with respect to each other, and produce 
a beat frequency when interfering on the detector. For 
a pair of particles situated at two different locations and 
moving with a velocity difference of Sv, the self-beat fre- 
quency induced by the Doppler shifts is q- Sv fA ITcI 
From here, it can be seen that measurements with q\\v 
will be influenced by these velocity-differences while mea- 
surements with q J_ v are not affected because the scalar 
products q- Sv are all zero. 

In the case of practical interest for our study, namely 
the situation in which the diffusion relaxation time is 
much shorter than the other relevant time scales, it is 
natural to assume that the scattered intensity fluctua- 
tions are well described by Gaussian statistics and that 
the intensity autocorrelation functions are related to the 
intermediate scattering functions by the Siegert relation- 
ship, 



g 2 (q,t) = l + (3-\gi{q,i)Y 



(6) 



where (3 is the speckle contrast which was, in the de- 
scribed experimental setup, of the order of 5%. In ref. 
[HI it is also mentioned that "it can be shown" that this 
relation is also valid for the case where the flow-induced 
deterministic motions dominate those due to diffusion, 
provided that the decay of the ISF is sufficiently short 
compared to the inverse local velocity gradient 7 _1 . This 
condition is almost always fulfilled with typical experi- 
mental parameters, and in particular is satisfied for our 
longitudinal flow scattering geometry (<f || v), which jus- 
tifies the use of the Siegert relationship in all the experi- 
ments performed here. 

In conclusion, the homodyne correlation functions 
measured in our XPCS experiments can be well described 
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by, 

l5i(<T,*)| 2 = \9i,nm\ 2 ■ \9i.r(t)\ 2 ■ \9i,sm\ 2 , (7) 

where the first factor (subscript d) is due to thermal dif- 
fusion, the second one (subscript t) to the transit time 
through the scattering volume, and the last one (sub- 
script s) is due to shear. Each of these factors are dis- 
cussed in the following. 



A. Thermal diffusion 

The intermediate scattering function for a non-flowing 
sample undergoing Brownian motion is a simple expo- 
nential decay with the relaxation rate determined by the 
diffusion coefficient D and scattering wavevector q, 



\9iA<i,t)\ 2 



exp 



-2DqH] 



(8) 



In this case, the dynamics is isotropic and does not de- 
pend on the direction of q, but only on its absolute value 
q. When a uniform shear rate 7 is applied, the diffu- 
sion is in general enhanced by the shear and becomes 
anisotropic, changing the formula to [14] : 



\9i,D(q,t)\ 2 = exp 




Here, q^ and q± are the components of q parallel and 
respectively perpendicular to the direction of the flow. 
In a transverse scattering geometry (<f _L v), <7||=0 and 
Eq. © reduces to Eq. ([5]). This is obviously not true 
when q\\ v, in which case the diffusion is enhanced by the 
shear. However, at the shear rates used here (see table[T]), 
the products tj are of the order of 10~ 3 and Eq. © 
provides only a minor correction to Eq. ©, which can 
still be considered to describe well the thermal diffusion 
of the colloids. 



B. Transit-time effects 

When the colloidal suspension flows through the beam, 
new particles enter the scattering volume, replacing par- 
ticles that leave on the other side. The "refilling" of a 
scattering volume of transverse size h, leads to a loss of 
correlation of the dynamic structure factor that can be 
characterized by a frequency ft T oc vo/h. 

The Deborah-number is a measure for the importance 
of this effect compared to the decorrelation due to ther- 
mal diffusion taking place on a timescale (Dq 2 ) . In 
this context, the Deborah number can be defined as, 



De = v (hDq 2 )' 



(10) 



Typical values for two values of q a re g iven in table |H 

It was demonstrated @, EE EE E3] that the contri- 
bution of the transit effect to the correlation function 
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FIG. 3: Correlation functions taken at no flow, Q — (sam- 
ple in a capillary) at all values of the scattering wavevector 
measured, q= 1.8, 2.1, 2.8, 2.5, and 5.3 xl(T 3 A -1 . The inset 
shows the normalized correlation functions shifted vertically, 
for clarity, by 0, 1, 2, 3, and respectively 4, and fits with ex- 
ponential decays - Eq. ©, solid lines. The main graph shows 
the correlation times r = 1 /Dq 2 (multiplied by q 2 in order to 
eliminate the q-dependence) as a function of q. 



is reflecting the beam profile, being "scanned" by the 
flowing particles. As the beam is defined by rectangu- 
lar slits, the incident intensity follows a ^ sm J r ) ^ dig. 

tribution, the squared Fourier transform of a rectangle- 
function. Due to the guard slits, which are aligned to 
suppress the higher order peaks, the beam profile at the 
sample position can be well approximated by a Gaussian 
form, leading to 



\gi,T(t)\ 2 oc exp [-(vtrt) 2 



Shear-induced effects 



(11) 



The shear-induced contribution to the correlation func- 
tion, |<?i sCSj t)\ 2 can be written as a sum over all pairs 
of particles or as a double integral over the scattering 
volume which can be approximated by a line of length 
2R, 



\9i, s m)f 
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(qt5v(ri,r 2 ))dridr 2 . 



The integral can be performed analytically for a uniform 
shear rate (Couette geometry) [13], leading to 



\9i,smf 



sin(gpof) 

q\\vot 



(13) 



In Appendix [A] Eq. (TT^J) was also solved analytically 
for a different flow geometry, a Poiseuille flow described 
by a parabolic velocity profile Eq. ([2]), which provides 
a better description for the data presented here. The 
resulting shear-induced response to the correlation func- 
tions can be written using the error function of a complex 
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argument, 

l<?i,s(<f,*)| 2 



Wq\\tvo 




AiqutvQ 



(14) 



From Eq. ([13]) and (|T4]) . a shear- induced frequency can 
be defined as 



v ■ q- 



(15) 



If the scattering alignment is supposed to be perfectly 
transversal, q_Lv, rs=0 and the XPCS correlation func- 
tions are unaffected by the shear-induced effects. 

The relative importance of the shear-induced effects 
compared to thermal diffusion is usually described by the 
Peclet number P — As it can be seen in tabled the 
Peclet numbers are very high even at low shear rates, and 
the shear-induced effects are the dominating contribution 
to the measured correlation functions if a non-transverse 
scattering geometry is used. 

Following the discussion above, it turns out that the 
relative influence of the shear can also be expressed in 
terms of a different dimensionless number which will be 
called here Shear number, given by the ratio between the 
characteristic diffusion time Td = 1/Dq 2 and the shear- 
time, t s = 1/T S = 1/q- v , 

q • v 



S = 



Dq 2 



(16) 



Values for the shear number S are given in table [J for 
a value of the longitudinal component of the scattering 
wavevector <7|| = 1.4-10~ 3 A -1 . The conclusion is the same 
as the one resulting from the evaluation of the Peclet 
numbers - the shear induced effects are the dominating 
ones if a non-transverse scattering geometry is used and 
it is hard or impossible to measure the thermal diffusion 
of the particles in such a scattering geometry. 

With the three components of the correlation func- 
tions determined by thermal diffusion, shear, and transit- 
time described by Eqs. ©, (Ti"4"|) , and respectively ([Til , 
the normalized intensity fluctuation autocorrelation func- 
tions measured in a XPCS experiment with a sample un- 
dergoing shear flow can be written as, 



lffi(q,*)l=G&(q,*)-l)//3 
exp 



-2D^(l-m^ *ll (7*)' 



r 



exp [-{v tT t) 2 ] 



16gi|iu 



liiq\\tv \ 



(17) 



As it will be shown in the following section, Eq. (TT7|) fits 
very well the measured correlation functions (Fig. 0]) 



IV. RESULTS &: DISCUSSION 

In the g-range probed here, the static structure fac- 
tor (FigfJl inset) of the colloidal suspension in glycerol 
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FIG. 4: Normalized correlation functions for a single wavevec- 
tor, q = 2.1 -10~ 3 A 1 and three different flow rates: a) Q = 2, 
b) Q = 10, and c) Q = 15 /il/h. The dashed line shows the 
fit of a simple exponential to the zero-flow (reference) corre- 
lation function for the same value of q. The solid lines show 
fits with (fl9|) and ((20)) as described in the text. 



exposed to a shear flow is isotropic and independent of 
the flow rate. This means that the flow does not change 
the spacial probability distribution of the particles. How- 
ever, the dynamic structure factor (intermediate scatter- 
ing function) is both flow-dependent and anisotropic. 

Fitting the data obtained from a non-flowing sample 
with a simple exponential decay - Eq. (|5|), as shown in 
Fig. [3] (inset), yields a mean diffusion coefficient of D = 
2.39 • 1CT 15 — . Using the Stokes-Einstein relation 



D 



k B T 
(rnrja ' 



(18) 



where k B is the Boltzmann constant and T = 293 K the 
temperature, the dynamic viscosity 77 can be calculated 
to be 0.839 Pa-s which corresponds to a left-over water 
content of 3% 

Using the obtained diffusion coefficient and viscosity, 
the flow can be completely characterized (see relevant 
values in tabled]). It can be concluded that the Reynolds- 
numbers are several orders of magnitude below the onset 
of turbulent flow. The Deborah-number - Eq. PU|) is q- 
dependent. At the smallest q- values measured, it is not 
possible any more to assume De <C 1 (even at low flow 
rates) and the transit-time effects must be considered 
for an accurate description of the correlation functions. 
At higher values of q, the Deborah numbers De become 
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TABLE I: Characterizing quantities of the flow of the sample. As a function of flow rate Q, the mean velocity vo, Reynolds- 
number Re, Deborah-number De (for two q- values), shear rate 7, Peclet number P, and Shear number (|16p are shown. 
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much smaller and the transit-time effects are less impor- 
tant. Increasing the beam size h offers another practical 
possibility to lower the Deborah number, if working with 
a lower speckle contrast is acceptable. The Peclet num- 
ber P and/or the q-dependent shear numbers S are much 
larger than unity even at the smallest flow rates, showing 
that the shear-induced effects are always important. 

The intensity fluctuation correlation functions were 
measured for several values of the scattering wavevec- 
tor q in both transverse (qJ-v) and longitudinal (<f|| v) 
scattering geometries and for several flow rates between 
and 80 /d/h and fitted using Eq. (TTT)) . For a transverse 
scattering geometry q\\ =0, and Eq. (TP?)) leads to 



92,±.{q,t) = l + /3-exp 



t 



exp [-K*) 2 ] , (19) 



with the diffusion time r = 1/Dq 2 and the transit- 
induced frequency vtr as fitting parameters (in addition 
to the speckle contrast (3 and the baseline value). The 
diffusion coefficient was obtained by fitting the low-flow 
transverse, q _L v correlation functions and was subse- 
quently fixed in the high-flow data which were used to 
determine the flow-dependent, transit-induced frequency 
24f Both parameters r and ist T were fixed to the values 
obtained from transverse scans when fitting the q\ \ v data 
with 



52,|| (g,*) = l + /3-exp 



-2l I , 



r . . 01 7T I erf (\/iTst) | 2 , . 

• exp [-K0 2 ] • ' 4r s / ' ( 2 °) 

from which a fitted shear-induced frequency T$ was ob- 
tained. 

Example fits for a single value of q and three different 
flow rates are shown in Fig. 2J The parameters resulting 
from the fitting procedure described above, vtr and Tg, 
are plotted in Fig. [5] as a function of the nominal flow 
velocity vq computed from the volume flow rate. The 
shear rate 7 was left in place in Eq. (T2TJ)) only for rea- 
sons of logical consistency with Eq. pTj) and ([9]), but its 
influence on the shape of the correlation functions is too 
small (for any reasonable values of the shear rate) to al- 
low any meaningful fitting. As it can be seen from Fig.|U 
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300 




FIG. 5: The fitted parameters versus the flow velocity of the 
sample; a) shows the results of fitting the transverse correla- 
tion functions with Eq. ([19"]) . Fits with Eq. ([20]) of the longi- 
tudinal correlation functions yields the values for the shear- 
induced frequency Ts displayed in b). 



the quality of the fits is quite remarkable, especially for 
q || v and for the higher flow rates, even if the only dom- 
inant parameter for these fits is Ts, the shear-induced 
frequency. 

From a practical point of view, in order to be able 
to measure the diffusive component of the particle dy- 
namics under flow, its relative importance must be high 
compared with that of the shear and transit-time induced 
effects. Keeping a transverse scattering geometry is the 
only method to limit the shear-induced effects, but in 
practice any small misalignment becomes the limiting 
factor at high-enough flow rates. If <j> is the angle between 
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FIG. 6: Dispersion relationships for the diffusion (Fd = Dq 2 ), 
shear (Fs = vo(j>q), and transit (y tr = vo/h) relaxation rates. 
The (example) values used here to estimate Fd, Fs, and v tr 
were D « 2.4- 10 5 A/s, ps 0.01 (0.5 deg), v = 1/im/s. The 
thick solid line shows the q region where the thermal diffusion 
dominates the dynamic signal and can be measured by XPCS. 
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FIG. 7: Correlation times (multiplied by q 2 ) as a function of 
q. At low flow rates the correlation times are unaffected by 
the transit time and/or the shear and the diffusive dynamics 
of the particles can be measured in a flowing solution. 



q and the direction orthogonal to Vp, the shear-induced 
relaxation rate - Eq. ijTij]) - becomes [lj|, 

T s = v Q <f>q. (21) 

To illustrate the impact of even small misalignments, 
the dispersion relationships for thermal diffusion Trj = 
Dq 2 , shear-induced effects, Eq. (f2"T]) . with <fi assumed to 
be 0.01 (corresponding to a misalignment of 0.5 deg) and 
vq=1 /im/s, are plotted in Fig. [6] together with the (q- 
independent) dispersion for the transit-time effects. In 
order to be able to measure the thermal diffusion of the 
particles, To must dominate over the other, flow-induced 
relaxation rates. For the values chosen in this example, 
this condition is fulfilled in the q-range for which the dif- 
fusion dispersion relationship is highlighted. From Fig. [B] 
it is clear that for any flow velocity larger than a few 
//m/s, the shear-induced relaxation rates Ts will shift up- 
wards which would "push" the possible q-range to higher 
values. Equivalently, working at a fixed value of q in order 
to obtain specific information about the dynamic prop- 
erties on a certain length-scale puts a condition on the 
maximum flow rate for which the correlation functions 
are still dominated by their diffusive components. At the 
same time, it is also clear that in the case of light-PCS, 
with smaller values of q, these conditions are harder to 
fulfill [3. 

From this analysis, it results that fitting the correla- 
tion functions obtained in transverse scans with Eq. (TIT))) 
is justified only at the lowest flow rates probed here (as- 
suming that it is difficult to achieve a vertical alignment 
better than 0=0.01). For such low flow rates, the cor- 
relation functions will not depend much on neither the 
finite transit-time nor on the shear time and fits with a 
simple exponential form will render the correct diffusion 
time and/or diffusion coefficient. 

The correlation times in Fig. [7] were obtained by per- 
forming fits to the transverse correlation functions with 
simple exponential forms. While at most of the higher 



flow rates probed here the influence of the flow on 
the correlation times is clear, at the low-enough flow 
rates (Q<2/il/h), the correlation times are still flow- 
independent and follow the rag 2 dispersion relation- 
ships expected for Brownian motion. 

Therefore, using XPCS on a flowing sample seems ap- 
plicable to measure the "intrinsic" dynamics of the sam- 
ple. While the scattering functions recorded with q || v 
are completely dominated by shear effects, these do not 
contribute (or contribute much less) to the spectra at 
qJ-V. 



V. SUMMARY 

The correlation functions from a colloidal suspen- 
sion undergoing shear flow were measured by homodyne 
XPCS. The data could be explained by computing the 
intermediate scattering function as a product of the scat- 
tering functions due to diffusion, transit, and shear. 

Our measurements suggest that it should be possible to 
measure samples under flow and obtain the diffusion co- 
efficient and/or particle size information, etc. For shear 
and Deborah- numbers < 0.1, the correction is negligi- 
ble and equation iJ5]) can be used. For higher Deborah 
and/or shear numbers, a more complicated convolution 
between different dynamic mechanisms takes place and 
the correlation functions are dependent on at least three 
important factors - thermal diffusion, transit-time, and 
shear. Measurements at shear and/or Deborah numbers 
above « 1 are completely dominated by shear-induced or 
(in some particular cases) transit time effects and not 
suitable any more for a precise extraction of diffusive 
properties. However, as both the De and S numbers 
are g-dependent, it is possible to tune the settings within 
some limits to access also higher flow rates. 

In practice, for a specific alignment a "calibration" pro- 
cedure (i.e. probing the dynamics of a colloidal suspen- 
sion as it is done here), should provide a graph such as 
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the one shown in Fig. [6] which provides a lower limit 
for q accessible for a given flow rate or, equivalently, the 
maximum flow rate that can be set while still being able 
to measure the diffusive dynamics at a specific value of q 
with the desired degree of accuracy. 

This method can therefore be used to measure the dy- 
namical properties of many soft and biological samples 
which suffer from beam damage, or to perform time- 
resolved studies in mixing devices. 
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APPENDIX A: SHEAR-INDUCED 
CORRELATION 



In the following, Eq. ([12]) is solved for a Poiseuillc flow, 
with a velocity profile described by and a closed- form 
expression for the shear-induced correlation \g% t s(q, t)\ 2 is 
calculated. 

The velocity difference between two particles located at 
different positions across the flow channel r\ and T2 can 

be written as bv(r\,T2) = 2 , the integrals over 



the two independent variables r% and T2 are equivalent 
and, after applying trigonometric addition identities, the 
double integral in equation (fl~2|) can be reduced to the 
integration over a single variable r, 



|0i,s(£*)| 2 = 7^ 



1 

4R2 



2 \ 



drcos ( ^q\\tv —j 



-i 2 



(Al) 



These integrals can be reduced to Fresnel integrals [20( , 
k,s(<7,i)| 2 - -r 1 — 

4q\\tv 




which are related by (using the imaginary number 



C(v)+iS(v) 
to the error-function, 

erf (u) 



e~* At. 



(A4) 



A normalized expression for the shear-induced correla- 
tion function |<7i,s(g, i)| 2 is, therefore given by 



\9i,s(q,t)\ 2 = 



16q\\tv 



or 



(A5) 
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